Measuring the Gluon Density in e + A Collisions: 
KLN CGC, DGLAP Glauber, or Neither? 
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We predict readily experimentally measurable differences in the diffractive cross section in the 
coherent exclusive photoproduction of J/ip mesons in e + A collisions at eRHIC and LHeC energies 
for nuclear gluon distributions assumed to 1) evolve in x with DGLAP dynamics and have a spatial 
distribution proportional to the Glauber nuclear thickness function and 2) evolve in x and b according 
to the KLN prescription of CGC dynamics. We find that CGC physics predicts that the nuclear 
gluon density widens significantly as a function of x yielding diffractive peaks and minima that 
evolve dramatically with x; on the other hand the DGLAP Glauber distribution yields peaks and 
minima constant in x. We also find that the dipole cross section at the level of two gluon exchange 
within the KLN parameterization of the CGC satisfies the black disk limit whereas this limit is 
violated when DGLAP evolution is used; the normalization of the diffractive cross section grows 
more slowly in x by several orders of magnitude when using the KLN parameterization as compared 
to the result when employing DGLAP evolution. 



I. INTRODUCTION 

Pcrturbativc quantum chromodynamics (pQCD) pre- 
dicts a nontrivial expansion in the size of the nuclear 
wavcfunction at small x due to the pcrturbativc power 
law tails of the gluon distribution near the edge of the 
nucleus [1—4]. Similarly, in order to not violate unitarity, 
the enormous growth in the gluon parton distribution 
function as x becomes small found via naive applica- 
tion of DGLAP evolution (see [5] and references therein) 
must be tamed by perturbatively-calculable saturation 
effects [2, 3]. However it is not yet clear from a theo- 
retical standpoint at what values of x these nontrivial 
changes in the dominant dynamics occur [3]. Addition- 
ally a quantitative theoretical understanding of exper- 
imental heavy ion data requires a quantitative under- 
standing of the initial geometry of a heavy ion collision. 
Certainly observables such as the azimuthal anisotropy 
of particles [6-8] is correlated with the anisotropy of the 
initial geometry; surprisingly the event-by-event fluctu- 
ations in the initial geometry also strongly affect these 
observables [9, 10]. In particular the viscosity to entropy 
ratio {rj/s) of the quark-gluon plasma (QGP) found by 
comparing hydrodynamics simulations to heavy ion col- 
lision data is directly related to the eccentricity of the 
initial thermal quark-gluon plamsa distribution that is 
evolved hydrodynamically. Currently the uncertainty in 
the initial thermal distribution due to the uncertainty 
in the importance of saturation effects in the initial nu- 
clear profiles is large enough that it is not clear whether 
the physics of the QGP is better described by leading 
order weakly-coupled perturbative quantum chromody- 



namics (LO pQCD) or by LO strongly-coupled anti-de- 
Sittcr/conformal field theory (AdS/CFT) methods [7]. 
An experimental measurement of the spatial gluon dis- 
tribution in a highly boosted nucleus, and hence the rel- 
evant physics in this kinematic range, would thus be a 
very interesting and important contribution to our un- 
derstanding of QCD. 
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FIG. 1: Leading order Feynman diagram for the exclusive 
vector meson production of a J/ip meson. 

Exclusive vector meson production (EVMP) in c + A 
collisions has been proposed as a channel for just such 
a measurement [12, 13]. In this Letter we will focus on 
the production of heavy vector mesons, in particular J/ip 
mesons. To leading order, EVMP of a J/ip meson oc- 
curs in an c + A collision when a photon emitted by 
the electron splits into a c-c pair which communicates 
with the gluon density in the highly boosted nucleus via 
a two gluon exchange and subsequently forms a J/ip me- 
son and nothing else (we will be interested here in co- 
herent EVMP, in which case the nucleus remains intact); 
see Fig. 1 for a visualization of the process. It is pre- 
cisely this two gluon exchange which yields a diffractive 
measurement of the gluon density in a nucleus; sec Fig. 2. 
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FIG. 2: The diffractive cross section da/dt from Eq. 1 for 1) the DGLAP evolved dipole cross section, Eq. 2, with gluons 
spatially distributed according to the Glauber thickness function of the Woods-Saxon distribution and 2) for the dipole cross 
section from the KLN model of the CGC, Eq. 8, for x = 10" 3 , 10" 4 , 10" 5 , and 10" 6 . The various a™ ax values shown explore 
some of the systematic theoretical uncertainty for the KLN CGC calculation, which is clearly much smaller than the difference 
between the results from DGLAP evolution and those from saturation physics. The yellow band represents the 1-a uncertainty 
(and the dashed black curve the central value) in the DGLAP Glauber results due to the uncertainty in the extracted LO 
MSTW gluon PDF [11]. 



Previous work [13] explored how modest changes in 
the Woods-Saxon distribution [14] of a nucleus might 
manifest themselves as changes in the diffractive peaks 
in EVMP if one assumes that the spatial distribution of 
gluons in a nucleus is proportional to the Glauber thick- 
ness function found from the Woods-Saxon distribution. 
That these modest changes do result in a visually obvious 
modification of the diffraction pattern motivated our fur- 
ther study, in which we consider whether two very differ- 
ent physical pictures of the gluon distribution in a highly 
boosted nucleus can be experimentally distinguished via 
EVMP: in particular we wish to compare the diffraction 
patterns that emerge when the gluon distribution 1) has 
normalization dictated by DGLAP evolution and spatial 
distribution given by the Glauber thickness function and 
2) is given by the KLN parameterization (see [15, 16] and 
references therein) of the Color Glass Condensate (CGC) 
(see, e.g., [3, 4] for a review). We choose to investigate 
these two ansatze of the gluon distribution in nuclei as 
they have been the dominant models used in heavy ion 
physics calculations to estimate the uncertainty in the 
viscosity to entropy ratio of the QGP produced at RHIC 
due to the uncertainty of the currently poorly constrained 
initial conditions in heavy ion collisions [6, 7]. 

It is worth taking a moment to comment on some 
common — yet confusing — terminology in the EVMP 
field. As mentioned above, to leading order the coher- 
ent production of a vector meson in an c + A collision 
involves a two-gluon exchange between the q-q pair and 
the nucleus. If one assumes that all two-gluon exchanges 
occur independently, then one may exponentiate the sin- 
gle two-gluon exchange result. Making this independence 
assumption is often referred to in the EVMP field as using 
"saturation" physics because the cross section is unita- 
rized via the exponentiation process. However this "sat- 



uration" does not refer to small-x evolution effects in the 
gluon distribution. For instance in the "IP-Sat" [17] and 
"b-Sat" [18] models, where "Sat" is short for saturation, 
the x evolution of the gluon PDF is effected through the 
use of the DGLAP equations. On the other hand, the 
"b-CGC" model [18] incorporates both the exponentia- 
tion of the two-gluon exchange and the CGC physics of 
the saturation of the gluon PDF. We note that, in prin- 
ciple, small-x evolution effects and exponentiation effects 
in the dipole cross section should become appreciable si- 
multaneously [19, 20]. In order to (hopefully) make the 
presentation more clear, and to simplify some of the nu- 
merics, we will not exponentiate the two-gluon exchange; 
we will present results using only the leading order two- 
gluon exchange in which the gluon PDF is given either 
via DGLAP evolution or from the CGC. Any subsequent 
reference to "saturation" in this paper will refer to the 
saturation of the gluon distribution function alone. 



II. FORMALISM 

Following [13, 17], the diffractive production of a vector 
meson from a photon scattering off a target is 
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where (V\j)t is the overlap of the vector meson wave- 
function and the transversely polarized virtual photon 
wavefunction — the contribution from the longitudinally 
polarized photon is zero as we are interested in Q 2 = 
photoproduction — and we used the photon-meson over- 
lap and Gauss-LC model for the J /ip wavefunction from 



[17] 1 , and A 2 = —t. da q q/d 2 b is the differential cross 
section for the interaction of the dipole with the target; 
its form depends on the physics assumptions we make 
for the nuclear gluon distribution, as we discuss in detail 
below. 



A. DGLAP Evolution in x, Glauber Distribution of 
Gluons in b 

If we assume that the two gluon exchange from the 
dipole to the nucleus occurs within an individual nucleon 
then 
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(2) 



where r is the size of the dipole, fi = ^(fio + C/r 2 ) is the 
relevant momentum scale for the dipole, xg is the gluon 
distribution function, and 



T(b) 
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(3) 



is the assumed spatial distribution of gluons in a nucleon. 
We use the MSTW parameterization of the gluon PDF 
[11]. As described in [21], /iq and C are free parameters; 
as in [17, 21, 22], we take fi = 1 GcV 2 and C = 4. 
From HERA data [23] the measured slope of dajdt yields 
B G « 4.25 GeV- 2 [17]. Then 
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where Jq is the usual Bcsscl function, 



T A (b) = j dz PA (y&+#) 
J d 2 bT A (b)=A, 



(4) 

(5) 
(6) 



is the usual thickness function, and p A is the density of 
the nucleus (here taken as the Woods-Saxon distribution 
of 197 Au with the usual R = 6.38 fm and a = 0.535 fm 
[24]) and 
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(7) 



1 Note that the normalization of the J/ip wavefunction in [17] is 
erroneously reported as a factor of 100 smaller than the correct 
value; one can readily see this by comparing with the normal- 
ization condition defined in [17] and with the results reported 
in [18]. It is surprising that this error was not noted in [18], in 
which the results found in [18] are explicitly compared to those 
in [17]. 
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FIG. 3: The two gluon dipole cross section, da q q/d?b from 
Eq. 2 for the DGLAP Glauber calculation and Eq. 8 for the 
KLN version of the CGC (multiplied by 50), which is propor- 
tional to the gluon density distribution in a nucleus probed 
by a q-q dipole of size r — 1 GeV -1 at x = 10 -5 . 



B. CGC Distribution of Gluons in x and b 

Alternatively we may view the nucleus as a whole and 
that the gluon distribution is found from the CGC. In 
this case 



da qq 7T 2 r 2\ I 2 
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(8) 



where xg A is the integrated gluon distribution function 
related to the unintegrated gluon distribution (UGD) <fi A 
by 



xg A (n 2 , Q 2 S ) 



d 2 kcj) A {k 2 



dk 2 cj) A (k 2 ,Q 2 s ) (9) 



The x and b dependence of the two-gluon exchange dipole 
scattering formula, Eq. 8, comes in implicitly through the 
x and b dependence of Q 2 [16], 



Q: 



2^2 

C F 



a s (Q 2 )xg(x, Q 2 s )T A (b), 



(10) 



where C F = (N 2 - 1)/2N C . 

In principle one determines the UGD via the JIMWLK 
evolution equations or, in the large- N c limit, the BK evo- 
lution equations (see [3, 4] and references therein). How- 
ever, instead of solving the full evolution equations many 
heavy ion physics calculations use instead the KLN pre- 
scription of the CGC (see, e.g., [15, 16]), which attempts 
to capture the main feature of CGC physics; in par- 
ticular, the KLN UGD becomes saturated at momenta 
on the scale of the saturation scale Q s . Because of its 
widespread use in heavy ion physics calculations and in 
order to simplify our own calculations we, too, will use 
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the KLN UGD, 
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(11) 

where k is an 0(1) parameter meant to represent higher 
order corrections to the UGD, and A = 0.2 GcV [16]. 

In principle k is set by comparing to known experi- 
mental observablcs such as the measured multiplicity at 
midrapidity at RHIC [25-27] or LHC [28, 29] or to the 
diffractive cross sections for protons measured at HERA 
[23]. However we found that the results from the lead- 
ing order multiplicity formula [15] are linearly dependent 
on the cutoff taken for a s , af ax . The KLN UGD it- 
self, though, is not nearly as sensitive to a™ ax , so the 
multiplicity prescription does not provide a robust way 
of setting k. We note in passing that the centrality de- 
pendence of the particles produced via the leading order 
CGC multiplicity formula using the KLN UGD's also de- 
pends on a™ ax . Perhaps the use of the next-to- leading 
order results in the UGD [30] and/or the production for- 
mula [31] will mitigate this dependence enough to make 
reasonable comparisons of CGC multiplicity to current 
data. Currently, though, there does not appear to be 
any quantitative estimate of the size of the dependence 
of the predicted CGC multiplicity as a function of cen- 
trality on ct™ ax . k also cannot be set by comparing to the 
proton diffractive cross section as the currently available 
data does not probe regions of x small enough such that 
Q 2 is a perturbative scale (at least when using the LO 
MSTW PDFs). In our calculations we will set k = 1. 

It is important to contrast the interaction of the dipole 
in the KLN CGC approach taken here, in which the q-q 
pair interacts with the entire nucleus, and the Glauber 
approach, in which the pair interacts with individual 
nucleons. By interacting with individual nucleons the 
diffractive cross section for the DGLAP Glauber model 
picks up an extra exponential suppression in t propor- 
tional to the square of the width of the nucleon, Be- 



lli. RESULTS 

In Fig. 3 we plot da q q/d 2 b 1 which is directly propor- 
tional to the gluon density probed by the heavy quark 
dipole, for Eq. 2 and Eq. 8, the DGLAP and KLN CGC 
distributions, respectively. The saturation physics of the 
CGC has resulted in a wider and flatter gluon distribu- 
tion than that from the Glauber treatment; the DGLAP 
growth of the small- a; gluon distribution — tamed by the 
saturation physics of the KLN CGC — leads to a signif- 
icant, nearly two orders of magnitude, enhancement in 
the cross section at a; = 10~ 5 compared to that found 
using the KLN CGC gluon distribution. It is worth not- 
ing that Fig. 3 shows that the KLN prescription for the 
CGC satisfies the black disk limit even at the level of 
two gluon exchange whereas application of LO DGLAP 



evolution leads to a violation of the black disk limit by 
an order of magnitude at x = 10~ 5 . 

We attempt to quantify the changes in both the nuclear 
gluonic width and density as a function of x in Fig. 4. In 
Fig. 4 (a) we show the quantity 61/2, which we define as 
the radius at which the dipole cross section reaches half 
its value at the origin: 



1 d(T q q 
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d 2 b 
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for the DGLAP and KLN CGC dipole cross sections. 
We note that even out to extremely small values of 
x ~ 10~ 13 , 61/2 from the KLN CGC continues to rise 
sublinearly with log(s); thus the implementation of the 
KLN CGC used here, with the MSTW gluon PDF, sat- 
isfies the Froissart bound [1]. Intriguingly this sublin- 
ear (as opposed to linear) growth in radius as a func- 
tion of log s is a surprise compared to other CGC pa- 
rameterizations [32]. In Fig. 4 (b) we show the depen- 
dence of the dipole cross section at b = on x for the 
DGLAP Glauber and KLN CGC models. Note the enor- 
mous growth of the dipole cross section as x decreases for 
the LO DGLAP-evolved gluonic density. This unitarity- 
violating enhancement is clearly reduced tremendously 
with the saturation physics of the KLN CGC. The yel- 
low band in the figure represents the 1-a uncertainty in 
the LO MSTW gluon PDF; the dashed black curve rep- 
resents the result using the central value of the LO gluon 
PDF [11]. 

In Fig. 2 we show the LO diffractive cross section 
for the EVMP of a J/ip in c + A collisions, Eq. 1, at 
x = 10~ 3 , 10~ 4 , 10~ 5 , and 10~ 6 when the gluon density 
grows in x and b according to DGLAP and Glauber over- 
lap or KLN CGC. As before, the yellow band describes 
the 1-a uncertainty in the LO MSTW gluon PDF, with 
the dashed black curve representing the central value. 
Several KLN CGC curves are plotted; they correspond 
to the results when the maximum cutoff for a s , a™ ax , is 
varied from 00 down to 0.5. Note that all previous fig- 
ures in this paper used OL^ ax = 00. For the various a™ ax 
curves in Fig. 2, the maximum value of the running cou- 
pling was set to a™ ax in: the dipole cross section, Eq. 8; 
the determination of the saturation scale, Eq. 10; and also 
in the KLN UGD, Eq. 11. While an interesting question, 
the influence of the uncertainty in the gluon PDF on the 
saturation scale is beyond the scope of this work. Clearly 
the KLN CGC diffractive cross section is not particularly 
sensitive to the specific ot^ ax chosen, which implies that 
higher order running coupling corrections to the result 
are small. The increase in the radial size of the gluon 
distribution as a function of a; shown in Fig. 4 (a) for the 
KLN CGC model manifests itself as a decrease in the 
spacings of the diffractive minima, At m j n j ma ~ l/b\/2, 
as one expects from a Fourier transform; on the other 
hand the positions in t of the maxima and minima of the 
diffractive cross section for the DGLAP Glauber dipole 
do not change as a function of x. 
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FIG. 4: (a) &1/2, the distance out from the center of the nucleus at which the dipole cross section is half its value at the center of 
the nucleus as a function of x. (b) da q q/d 2 b\b=o, the value of the dipole cross section evaluated at the center of the nucleus, as 
a function of x. The yellow band represents the 1-a uncertainty (and the dashed black curve the central value) in the DGLAP 
Glauber results due to the uncertainty in the extracted LO MSTW gluon PDF [11]. Both (a) and (b) are evaluated at dipole 
size r — 1 GeV -1 . 



The drastically faster increase in the gluon density 
from the DGLAP evolved PDF seen in Fig. 4 (b) results 
in a cross section that increases much faster as a func- 
tion of x than for the KLN CGC case. As was shown in 
[22] 2 the incoherent cross section, in which the nucleus 
breaks up, begins to dominate the total diffractive cross 
section by t ~ 0.02 GeV -2 . It is likely that the t depen- 
dence of the incoherent EVMP of the two models will 
be different, although we do not provide a quantitative 
estimate here: the decrease in cross section as a function 
of t for the DGLAP Glauber model will be enhanced by 
exp(— Bq t) due to the assumption that the heavy quark 
dipole interacts with individual nucleons. And in the case 
of coherent scattering shown in Fig. 4, one can discern 
a stronger t dependence in the DGLAP Glauber results 
due precisely to the extra exp(— Bq t) factor that results 
from treating the nucleus as a collection of individual 
nucleons. More importantly, the much larger gluon den- 
sity yields a particularly noticeable difference at t = 0, 
where possible nuclear breakup effects are negligible: the 
DGLAP Glauber case is an order of magnitude larger 
than the KLN CGC case at x = 10~ 3 and is a full two 
orders of magnitude larger at x = 10~ 6 . Even with 
the very large PDF uncertainties as x decreases, there is 
a clear increase in the coherent diffractive cross section 
for the DGLAP Glauber dipole compared to the KLN 
CGC dipole. Note that the enormous normalization dif- 
ferences seen in Fig. 4 (b) between the DGLAP Glauber 



2 Figure 8 in [13] also shows that the incoherent process quickly 
dominates the coherent one as a function of t, although we note 
that there was an error in the calculation of the figure and that 
the curves plotted do not correspond to the equations in the text 
of the paper. 



and KLN CGC dipoles for the most likely dipole size of 
r = 1 GeV -1 for the photon- vector meson overlap do 
not directly translate into as large normalization differ- 
ences in da/dt due to the integration over all dipole sizes, 
r. 



IV. CONCLUSIONS AND DISCUSSION 

An enormous wealth of information on the gluonic 
structure of highly relativistic nuclei can be found us- 
ing exclusive vector meson production. In particular we 
investigated the experimental signatures of the coherent 
scattering of a cc dipole onto a nucleus that results in an 
intact nucleus and a J/ip meson in e + A collisions at 
cRHIC and LHeC energies. We found that the diffrac- 
tive cross section will readily experimentally differentiate 
between the two common initial highly boosted nucleus 
prescriptions used in heavy ion physics phenomenology: 
1) the gluon density is found using DGLAP evolution 
and its spatial distribution is assumed to be proportional 
to the at-rest Glauber nuclear thickness function and 2) 
the gluon density and distribution is given by the KLN 
parameterization of the CGC. In particular there is the 
exciting possibility of literally watching a nucleus grow 
with center of mass energy as the positions in t of the 
minima and maxima in the diffractive cross section for 
the saturation physics calculation depend quite strongly 
on log(x). On the other hand the DGLAP Glauber model 
yields a nucleus of constant size as a function of 2; the 
positions in t of the diffractive minima and maxima do 
not change as a function of x. At the same time one is 
determining the width of a nucleus in e + A collisions, 
one will also measure the x dependence of the normaliza- 
tion of da/dt. Due to the explosion of small-a; gluons the 
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DGLAP Glauber approach yields a normalization that 
rapidly increases as a function of x; additionally the t 
dependence of the DGLAP Glauber da/dt is also quite 
strong as it is proportional to exp(— Bq t) due to the as- 
sumption that the q-q dipole interacts with individual 
nucleons. Conversely the KLN CGC dipole description 
does not have a strong x dependence in its normaliza- 
tion due to its inclusion of saturation effects; similarly, 
the interaction of the dipole with the whole nuclear glu- 
onic wavefunction yields a weaker t dependence than is 
displayed by the DGLAP Glauber results. 

It is clear that, at the very least, the striking difference 
between the x dependence of the peaks and minima from 
the DGLAP Glauber model and the KLN CGC model are 
robust: these differences will persist should we use even 
more sophisticated models of these two physical pictures; 
the x dependence of the peaks and minima will persist 
should we attempt to approximate multiple scattering 
within the nucleus by exponentiating the dipole cross sec- 
tion, should we use a less approximate CGC calculation 
such as is found in [30] , or should we examine the results 
from other vector mesons such as the <f> or p. We re- 
grettably leave the quantification of the diffractive cross 



section for these more sophisticated physical models and 
additional vector mesons for future work. Exponentiat- 
ing the two-gluon exchange cross section will reduce the 
enormous growth in the diffractive cross section in the 
DGLAP Glauber picture compared to the CGC case; we 
suspect this reduction will not be too large, although we 
also leave the quantification of this reduction to future 
work. 
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